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We propose that a half semiconducting state can exist in trigonal-prismatic transition metal 
dichalcogenide (TMDC) monolayers of d^ configuration. In that state both electrons and holes are 
spin polarized and share the same spin channel. On the basis of hybrid density functional theory, 
we predict in particular that VS2 monolayers are half semiconductors with a direct band gap. 

Moreover, we find that the conduction electron spin orientation of VS2 switches under moderate 
strain. Our predictions thus open up intriguing possibilities for applications of VS2 in spintronics 
and optoelectronics. Our analysis of trigonal-prismatic group-V MX2 (M=V, Nb, Ta; X=S, Se, Te) 
monolayers reveals a broad diversity of electronic states that can be understood qualitatively in 
terms of localization of d electrons. 
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Two-dimensional transition metal dichalcogenides 
(TMDCs) have been the subject of intense research over 
the past decade because of their unusual electronic prop¬ 
erties and potential applications in future nanoelectronic 
devices [ll-Q- The prototype of such materials is M0S2, 
which undergoes a transition from an indirect band gap 
in the bulk to a direct band gap in the monolayer [^. 
The direct-band-gap transition, together with its valley 
and spin selective properties, makes the material very 
promising for new generation electronics d, @ . Related 
properties were also observed in other groim-VI dichalco¬ 


genides such as MoSe2, WS2 and WSe2 [THSl- In compar¬ 
ison, the central topic in group-V TMDCs is electronic 
instabilities [IMl- Owing to complicated electron- 
electron and electron-phonon interactions, this series of 
compounds displays a rich phase diagram which includes 
metal, charge density wave (CDW), Mott insulator, and 
even superconductor [l^ . In addition, their electronic 
structure appears to be highly sensitive to temperature 
and pressure, which can induce CDW or metal-insulator 
transitions [ij. The complexity of the electronic struc¬ 
ture of these materials holds promise of more surprising 
physics yet to be discovered. 


In this Letter, we study the electronic structures of 
trigonal-prismatic MX2 (M=V, Nb, Ta; X=S, Se, Te) 
monolayers based on first-principle calculations at the 
level of the hybrid functional, and propose that the 
ground state of VS2 monolayers is a correlation-driven 
half semiconducting state. In this intriguing half semi¬ 
conducting state — examples of which are rare — both 
the valence band maximum (VBM) and the conduction 
band minimum (CBM) are spin-polarized, and their spin 
orientations are identical. Moreover, the spin sign of the 
CBM in VS2 is found to switch under moderate strain, 
undergoing a phase transition from a half semiconductor 
to a magnetic semiconductor. This opens up intrigu¬ 
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FIG. 1 : (Color online) Energy band diagrams for trigonal- 
prismatic group-V MX2 monolayers with increasing spin split¬ 
ting of valence band. Depending on U/t, four different mean- 
field phases are expected: (a) nonmagnetic metal, (b) mag¬ 
netic metal, (c) magnetic semiconductor, and (d) half semi¬ 
conductor. 


ing possibilities for applications in novel spin based elec¬ 
tronics. Our analysis of other Group-V MX2 monolayers 
(M=V, Nb, Ta; X=S, Se, Te) reveals a broad range of 
correlation-driven phases such as ferromagnetic metallic, 
ferromagnetic semiconducting states, and half semicon¬ 
ducting states (Fig. [T]). We demonstrate that this rich 
phase diagram — interesting both from a fundamental 
as well as an applied perspective — has an intuitive in¬ 
terpretation in terms of localization of d electrons. 

A trigonal-prismatic MX2 monolayer is composed of 
one triangular transition-metal sublattice sandwiched be¬ 
tween two sheets of chalcogen atoms. Chalcogen atoms 
are below and above the centers of metal triangles, re¬ 
sulting in a six-fold trigonal-prismatic coordination for 
the metal , cf. Fig. [ 2 j The coordination splits the 
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TABLE I: Electronic states for trigonal-prismatic MX2 
(M=V, Nb, Ta; X=S, Se,Te) computed at the level of hy¬ 
brid fnnctional DFT. * Although VSe2 and VTe2 are magnetic 
semiconductors, their VBM and CBM arise from splitting of 
different bands, similar to half semiconductors, unlike Fig. [T]:. 
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FIG. 2 : (Color online) Top view and side view of trigonal- 
prismatic MX2 monolayers. (M atoms in red and X atoms in 
yellow). 


metal d orbitals into A[ E' {dxy,dx^-y^), and E" 

[dxz,dyz) manifolds [l^. Governed by the overlap with 
chalcogen orbitals, their mean energies lie in ascending 
order of E\ and E". The d^2, dxy, and dx^-y'^ or¬ 
bitals further hybridize with each other owing to common 
reflection symmetry, contributing a narrow valence sub¬ 
band and a wider conduction subband [Sin. In other 
words, the low energy bands of MX2 are dominated by 
the onsite hybrid and the inter-site overlap of the dz2 , 
dxy, and dx;2_y2 orbitals of the 2 D metal sublattice. A 
tight-binding model taking into account these three or¬ 
bitals already captures well the dispersion of the low en¬ 
ergy bands — especially the valence band 18 |. Also re¬ 
cent experimental investigations provide evidence for the 
metal-metal interactions to define the electronic struc¬ 
ture of such 2 D systems For these reasons, a multi¬ 
orbital Hubbard model on the triangular lattice is a good 
candidate for describing the strongly correlated physics of 
d^ MX2 monolayers. However, to gain a qualitative un¬ 
derstanding of the emergence of the various correlation- 
driven phases shown in the schematic diagrams of Fig. [TJ 
we use a simpler, single band Hubbard model: 


H = -tY^ {ct^Cja + c^aCia) + U ( 1 ) 

i,j,cr i 

Here, t is the nearest-neighbor hopping integral, and 
U is the on-site coulomb potential. cf^{cia) creates (an¬ 
nihilates) an electron of spin a at site i, and riia is the 
number operator for spin a at site i. U and t are two 
competing parameters that parametrize the localization 
and the mobility of the strongly correlated electrons, re¬ 
spectively. The onsite repulsion U is responsible for the 
spin splitting of the narrow d band. As illustrated in 
FiglU on a mean-field level, there are four possible band 
structure configurations, depending on the magnitude of 
the spin splitting. Without any spin splitting (Fig. [ 1 ] 
a) the system is a half-filled nonmagnetic metal. In the 
second case (Fig. [Ud) the two spin components of the 


narrow d band are partially separated and the system is 
a magnetic metal. As the spin splitting starts exceed¬ 
ing the bandwidth, such that the two spin components 
are fully separated (Fig. (T}:), a magnetic semiconductor 
results. Finally, when the spin splitting of the narrow 
d band becomes larger even than the crystal field split¬ 
ting, so that one spin component of the narrow d band is 
above the CBM, which is typically about 1.0 eV higher 
than the VBM for TMDCs, the system is a half semi¬ 
conductor (Fig. Id): both the VBM and the CBM are 
spin-polarized and have the same spin orientation. The 
first three states are frequently reported while the fourth 
one is rare in nature. Interestingly, as our calculations 
show, d^ TMDCs span the entire gamut of phases shown 
in Fig. [U including the half semiconducting state, which 
we predict to exist in VS2. We note that a large spin 
splitting comparable to the crystal field splitting does 
not always result in a half semiconductor: the VBM and 
CBM originating from two different spin-split bands may 
end up having different spins (as in the case of VSe2 and 
VTe2; see below). 

We make the above predictions for the electronic struc¬ 
tures of trigonal-prismatic MX2 (M=V, Nb, Ta; X=S, 
Se, Te) monolayers based on first-principles calculations. 
(We note that although for VSe2 and VTe2 the trigonal- 
prismatic configuration is not the ground state, we in¬ 
vestigate it here for its theoretical interest.) The cal¬ 
culations are performed using density functional theory 
as implemented in the Vienna ab initio simulation pack¬ 
age (VASP) 0 - The exchange-correlation functional 
is considered at the level of the generalized gradient 
approximation (GCA) [2^. The screened hybrid func¬ 
tional method proposed by Heyd, Scuseria, and Ernz- 
erhof (HSE) [ 21 j is used. The single-particle equations 
are solved using the projector-augmented wave (PAW) 
method [H, with a plane-wave basis and a cutoff en¬ 
ergy of 500 eV. All the calculations are carried out in a 
slap model with the primary surface cell without consid¬ 
ering any CDW distortion, since we focus on the electron- 
electron interaction only. Our results are summarized in 
Table HI While all three VX2 (X=S, Se, and Te) mate¬ 
rials show large spin splitting comparable to the crystal 
field splitting, only VS2 has the same spin orientation 
for both CBM and VBM, and, therefore, is a half semi- 
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FIG. 3 : (Color online) Spin-resolved (a) band structures and 
(b) density of states for VS2 computed with hybrid DFT. The 
presented V 3 d density of states are contributed by d^2, d^y, 
and d^2_y2 orbitals. 



conductor. In trigonal-prismatic VSe 2 and VTe 2 they 
have different spin orientations (see Supplementary ma¬ 
terials). Like NbSe 2 , NbTe 2 , TaSe 2 , and TaTe 2 they are 
thus ferromagnetic semiconductors, while NbS 2 and TaS 2 
are ferromagnetic metals. 


Half semiconducting state and direct band gap in ¥ 82 - 
According to our total energy calculations, the VS 2 
monolayer is in a ferromagnetic ground state, consistent 
with both prior calculations 2 ^ and experimental obser¬ 


vations [2^, [ 23 . The ferromagnetic state is more stable 
than both the nonmagnetic and the antiferromagnetic 
states (spins of nearest neighbors are opposite), by 367 
meV and 689 meV, per cell, respectively. The magnetic 
moment per unit cell is calculated to be 1.0 piB. Each V 
atom contributes 1.30 pB, while each S atom contributes 
-0.15 plB. As is seen from Fig. O the electrons near both 
the VBM and the CBM have the same spin orientation, 
making VS 2 monolayers moreover half semiconductors. 
The highest valence subband is made up of spin-up hy¬ 
brid states of dz^, d^y, and dx 2 _y 2 . The Coulomb in¬ 
teraction causes a ^^2.5 eV spin splitting, pushing the 
down-spin counter part of this band into the conduction 
band. The CBM consists of spin-up electronic states of 
the higher hybrid band. As a result, the total density of 
states matches well the half semiconducting state illus¬ 
trated by the schematic Fig. [TJl. In addition, both the 
VBM and CBM are located at the K point, with a direct 
band gap of ^1.1 eV. The highest occupied state at K is 
93 meV higher than the one at the T point. The lowest 
unoccupied state at K is 316 meV lower than that of the 
opposite spin at a local minimum near M. The direct 
band gap of the VS 2 monolayer is similar in magnitude 
to that of the isostructure M 0 S 2 monolayer, and it is well 
in the region of the visible spectrum — of great interest 


FIG. 4 : (Color online) Spin-resolved electronic structures for 
(a) 3 % contracted and (b) 3 % expanded VS2. 


for optoelectronic applications. We arrive at similar con¬ 
clusions by performing additional, GGA-I-U calculations 
(see Supplementary materials). 

Experimentally, VS 2 monolayers have not yet been re¬ 
alized. So far only thin films of VS 2 have been reported 
[ 2 ^ . They are metallic, as a result of being in the oc¬ 
tahedral phase and because of the interlayer interaction 
[ 2 ^. Further experimental efforts are needed to confirm 
the intriguing properties we predict for VS 2 monolayers. 

Strain induced phase transition in VS 2 monolayer. Me¬ 
chanical strain provides a means to tune the electronic 
structure of 2D electronic systems by controlling the 
spacing between atoms. We, therefore, next examine the 
effect of strain on the electronic structure of VS 2 . Ac¬ 
cording to our calculations, compressive strain tends to 
decrease the energy of the half semiconducting state and 
destroy it while tensile strain tends to stabilize it. Fig. 
m shows the band structure and density of states of a 
VS 2 monolayer under 3% strain, i.e. the lattice constant 
was contracted or expanded by 3%. Under compressive 
strain, the atomic spacing of V atoms decreases. This 
increases the orbital overlap and hence the hopping inte¬ 
grals t of the d electrons. As a result, the width of the 
narrow d band is increased owing to the electron delo¬ 
calization and the spin splitting is reduced. The system 
transitions into the normal ferromagnetic semiconduct¬ 
ing state due to the lowering of the energy of the down 
spin counterpart of the valence d band (Fig. [4^). The 
lowest unoccupied state then is a down spin state near 
M, which has energy 327 meV lower than the lowest 
state at iT, which is spin-up. According to our calcula¬ 
tions, the critical strain for this transition is 2%. This is 
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a quite moderate, experimentally achievable value. On 
the other hand, tensile strain enlarges the atomic spac¬ 
ing and enhances the localization of the metal d elec¬ 
trons at the lattice sites. That localization in turn leads 
to a larger onsite interaction, which enlarges the energy 
splitting between opposite spins. Consequently, the half 
semiconducting state is enhanced as the down spin coun¬ 
terpart of the valence subband is lifted deeper into the 
conduction band (Fig. IDd). Our results reveal that the 
spin splitting of the d electronic states in VS2 is sensitive 
to the atomic spacing of V atoms owing to electron local¬ 
ization. Thus, the spin sign of the conduction band can 
be switched by moderate, compressive strain. This phe¬ 
nomenon has direct application in spin-based electronics. 

Electronic phases of group- V MX2 monolayers. While 
VSe2 and VTe2 have similar spin splitting as VS2, push¬ 
ing the down spin counterpart of the d valence subband 
into the conduction band, the VBM in those cases is not 
located in the spin-up d subband, but in a spin-split sub¬ 
band arising from deeper p bands. The VBM and CBM 
thus have different spin signs (see Supplementary mate¬ 
rials), such that these materials are unconventional fer¬ 
romagnetic semiconductors. On the other hand, NbS2 
and TaS2 are ferromagnetic metals, and NbSe2, NbTe2, 
TaSe2, and TaTe2 are normal ferromagnetic semiconduc¬ 
tors. This series spans the entire range of strongly cor¬ 
related phases illustrated in Fig IT] A summary of the 
phases is given in Table H) 

There is a trend for small metal and large chalcogen 
components to stabilize the phase on the large U/t side 
of our Hubbard mean-field schematic Fig. [l] which can 
be attributed to electron localization. Consider MX2 of 
the same metal but different chalcogen components X. 
In this case, the orbital properties of the metal compo¬ 
nents are fixed parameters, while the hopping integral is 
determined by the distance between the transition metal 
atoms. As shown in Fig. [ 5 ^, the ratio of metal diame¬ 
ter to the cell constant, c?M/a, decreases with increasing 
chalcogen radius from S to Te. Since this ratio quali¬ 
tatively measures the overlap of the d orbitals, this re¬ 
sults in a narrowing of the d band and an increase of 
the onsite interaction, i.e. a larger spin splitting. The 
phase variation and the trend of the energy gain by fer¬ 
romagnetic order, AE = EpM — Enm (Fig. [SJd) for 
MX2 with different chalcogen components can be under¬ 
stood by these simple considerations. If one considers 
instead different metal components keeping the chalco¬ 
gen component fixed, one can draw a similar conclusion 
(see Supplementary materials). 

The above results provide strong evidence for the 
strongly correlated nature of electronic states in group- 
V TMDCs. So it is not surprising that contradicting 
predictions exist in the literature. Based on the GGA 
method. Pan recently reported that yX? monolayers are 
all semimetallic and ferromagnetic [ 29 |. Wasey et al. 
predict VS2 to be a semimetal, and VSe2 and VTe2 to 



FIG. 5 : (Color online) (a) The ratio of M diameter to the cell 
constant, dAf/a, and (b) the difference between the energies 
Efm of the ferromagnetic and Enm of the nonmagnetic states 
as a function of X radius for MX2. (M=V, Nb, Ta; X=S, Se, 
Te.). 


be semiconductors with indirect band gaps 


Li et 


al. also predicted a semiconducting state for VSe2 [ 31 1 . 
These variations in the results can be addressed by asking 
the following question: are the opposite spin components 
of the narrow d band in contact or fully separated? It is 
hard to give a conclusive answer at the level of the GGA 
method since the calculated band structure is close to the 
transition between the two. Our calculations based on 
the hybrid functional clearly demonstrate the existence 
of the half semiconducting state. 

It should be noted that the predicted correlation driven 
ferromagnetic and insulating states are hard to observe in 
experiments as they are easily corrupted by doping or dis¬ 
order (sil - fs^l . The ferromagnetic, semiconducting state 
is based on a correlation-driven electronic configuration 
where the down spin component of the highest valence 
subband is a few eV higher than its up-spin counterpart 
and unoccupied. The energy gain over the nonmagnetic 
and metallic states from the spin exchange, however, is 
only a hundred meV or smaller. Occupation of the down- 
spin band would thus be energetically unfavorable, and 
doping collapses the system into the nonmagnetic and 
metallic state. Especially for NbX2 and TaX2 the stabi¬ 
lization energy for ferromagnetic order and the semicon¬ 
ducting state is on the order of only tens of meV, making 
these states likely too fragile to be observed. However, 
it is possible for the ferromagnetic and semiconducting 
states to survive in VX2 since they have relatively large 
stabilization energies on the order of a hundred meV. 
This may explain why ferromagnetism has been observed 
in VS2 multilayers. 


In summary, we predict that the electronic states of 
d^ TMDG monolayers span a rich variety of correlation- 
driven phases including magnetic metal, magnetic semi¬ 
conductor, and half semiconductor. Our predictions are 
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based on hybrid functional DFT calculations. Most im¬ 
portantly, our calculations demonstrate that the VS 2 
monolayer is a direct band gap half semiconductor, the 
electron spin orientation of which can be switched by 
moderate strain. Such easily tunable phase transitions 
are highly desirable and hold great promise for appli¬ 
cations in spintronics. Our study moreover reveals the 
importance of the role of electron localization in the elec¬ 
tronic structure of d^ TMDCs, which is of fundamental 
interest in the theoretical understanding of these materi¬ 
als. We anticipate our findings to inspire further experi¬ 
mental and theoretical work, particularly with a view to 
application in novel electronics. 
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Note added: Shortly before completion of this 
manuscript we became aware of complementary work on 
VS 2 , which confirms that the trigonal prismatic phase of 
the material studied here is indeed the ground state at 
the magnitude of strain considered here [35 1. 
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SUPPLEMENTARY MATERIALS 

GGA + U calculations for US 2 . In addition to hybrid 
calculations, we also carried out GGAj-U calculations 
to examine the effect of Hubbard U on the electronic 
structure of the VS 2 monolayer. According to our band 
structure calculations, the spin splitting of the narrow 
d band gradually when the onsite Goulomb potential U 
is increased from 0 to 4.0 eV (Fig. [5]). The energy of 
the down spin component of the valence subband con¬ 
tinually increases from that of its spin up counterpart 
until it is well above the CBM. As a result, the system 
undergoes a series of transitions, from semimetal to fer¬ 
romagnetic semiconductor to half semiconductor, as pro¬ 
posed in Fig. [TJ For a certain value of U (^3.5 eV), 
the GGA-I-U method results in the same prediction of 
a half semiconducting state as made by the HSE func¬ 
tional. The HSE functional includes some short range 
Hartree-Eock exchange, which partially cancels the self¬ 
interaction emerging in GGA and results in an increased 
localization of d electrons. The localization of electrons 
on atomic sites thus increases the onsite interaction. On 
the other hand, the GGA-I-U method imposes an extra 
onsite interaction, leading to larger spin splitting. The 
consistency of the two approaches demonstrates the criti¬ 
cal role of correlation due to localization in the electronic 
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FIG. 6: The low energy band structure of VS 2 as a function 
of U between 0 and 4.0 eV (0.5 eV per step). The thickness 
of the line indicates the value of U. 



FIG. 7: Energy difference between the ferromagnetic and non¬ 
magnetic states of monolayer VS 2 as a function of strain. 

structure of MX 2 . 

Strain and ferromagnetic ordering. In Fig. [71 we show 
the energy gain of ferromagnetic order versus the non¬ 
magnetic states EpM — E^m as a function of in-plane 
strain for monolayer VS 2 . As the strain varies from -4% 
(compressive) to 4% (tensile), the energy gain increases 
almost linearly from 0.30 eV to 0.41 eV. This demon¬ 
strates explicitly that compressive strain decreases the 
energy gain of ferromagnetic order and tends to destroy 
the half semiconducting state. On the other hand, ten¬ 
sile strain tends to enhance ferromagnetic ordering and 
hence the half semiconducting state. 



'^FMK r MK r 


FIG. 8: Spin-resolved band structures for monolayer VSe 2 
and VTe 2 . 


Unconventional magnetic semiconductors, VSe 2 and 
VTe^. In Fig. [51 we show the spin-resolved band struc¬ 
tures for trigonal-prismatic , monolayer VSe 2 and VTe 2 , 
respectively. The down spin counter part of the valence d 
band is above the CBM in these two compounds. How¬ 
ever, in both cases, the VBM is not located in the up 
spin d band, but in a spin split subband arising from 
deeper p bands. Therefore, VSe 2 and VTe 2 are unconven¬ 
tional magnetic semiconductors which differ from both 
half semiconductors like VS 2 and magnetic semiconduc¬ 
tors like NbSe 2 , NbTe 2 , TaSe 2 , and TaTe 2 . For VSe 2 , the 
VBM, made up of p bands, is only 40 meV higher than 
the highest d states at K. It may be possible to tune it 
into a half semiconducting state. However, for VTe 2 , the 
VBM is 731 meV higher, deviating greatly from a half 
semiconducting state. 

Effect of metal atoms on localization and ferromagnetic 
ordering. The phases in Table HI can be understood in 
terms of electronic localization on the metal atoms. In 
the main text we have shown how the different phases for 
the same metal atom with different chalcogen atoms can 
be understood within this picture. Here we consider the 
phase differences among MX 2 materials with the same 
chalcogen atom, but different metal atoms. The diame¬ 
ters of the group-V transition elements, V (3d), Nb (4d), 
and Ta (5d), are 3.42, 3.96, and 4.00 A, respectively. 
According to our calculation, the cell constants for their 
compounds with sodium, for example, are 3.17, 3.35, and 
3.35 A respectively. As a result, the dpi/a. ratio for VS 2 
(1.08) is much smaller than for NbS 2 (1.19) and TaS 2 
(1.20). Similar considerations apply for monolayers with 
Se or Te components. These comparisons show that the d 
electrons in VX 2 are more localized than those in NbX 2 
and TaX 2 . Consequently, all VX 2 monolayers are well 
stabilized in the semiconducting state, while NbX 2 and 
TaX 2 fall into the critical region of metal to insulator 
transition. Also, the stabilization energy of ferromag- 
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netic order for VS 2 , A.E, is 367 meV, much larger than than that in NbS 2 and TaS 2 owing to stronger electronic 
those of NbS 2 (122 meV) and TaS 2 (64 meV), and it fol- localization, 
lows the same trend: the spin splitting in VS 2 is larger 


